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mally crossed diagrams the pair-pair correlations are considered. It is shown 
that the recently proposed asymmetric Bethe-Salpeter equation to avoid un- 
physical repeated collisions is derived as a result of the hierarchical dependen- 
cies of correlations. Exceeding the parquet approximation we show that an 
asymmetry appears in the selfconsistent propagators. This form is superior 
over the symmetric selfconsistent one since it provides the Nambu-Gorkov 
equations and gap equation for fermions and the Beliaev equations for bosons 
while from the symmetric form no gap equation results. The selfenergy di- 
agrams which account for the subtraction of unphysical repeated collisions 
are derived from the pair-pair correlation in the three-particle Green's func- 
tion. It is suggested to distinguish between two types of selfconsistency, the 
channel-dressed propagators and the completely dressed propagators, with 
the help of which the asymmetric expansion completes the Ward identity 
and is (^-derivable. 
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1 Introduction 

Though pairing and condensation phenomena like superconductivity belong 
to the most exciting and long studied effects, their many-body theoretical 
treatments lacks often a systematic justifications. For example a particle- 
number non-conserving assumption is often accepted by anomalous propaga- 
tors in order to obtain the desired structure of Green's functions. Systematic 
many-body approximations like the ladder T-matrix [l][2l[3j|4l[5] leading to 
the Bethe-Salpeter equation had to face the puzzle that the selfconsistent 
propagators describing the many-body medium have to be assumed in asym- 
metric form in order to obtain the gap equation. This was first observed by 
Kadanoff and Martin [1] and this approximation has been used mainly in the 
theory of superconductivity [3J[5JE] . For an overview about different T-matrix 
approximations see [7]. Though any T-matrix based on the Bethe-Salpeter 
equation becomes unstable at the critical temperature [51 where the T-matrix 
diverges, the corresponding selfenergy fails to describe the superconducting 
gap [S] if we do not assume an asymmetry in the internal propagators. 

Recently a way has been proposed to obtain such asymmetry by demand- 
ing that particles should not interact with the same state repeatedly |2| . Such 
non-physical repetition leads to an incorrect description of the coherent part 
of the scattering [101 , and an incorrect single-particle Green function G. These 
non-physical repeated collisions destroy the superconducting gap in the en- 
ergy spectrum of G. Using the Soven scheme [IT the repeated collisions can 
be avoided [Qj and the gap in the superconducting state can be described as 
well as the normal state. This subtraction scheme can be written in a form 
equivalent to the anomalous propagator approach such that these propaga- 
tors appear as consequence of the theory and need not to be assumed ad-hoc 
|12) . Though in this way a physical understanding of a deficiency in selfcon- 
sistent many-body expansion has already been achieved, it is desirable to see 
how a systematic expansion of correlations can account for such corrections. 

The necessity to introduce asymmetric propagators can be seen from an- 
other physical system of dense correlated plasmas. There the motion of ions 
are lowered by the dynamics of the screening cloud since a deceleration force 
appears which is caused by the deformed screening cloud surrounding the 
charge. This effect known as Debye-Onsager relaxation effect has been de- 
voted much interest [T51[H[ [15.16.17.18,19,20 which was first derived within 
the theory of electrolytes [?lu22n23u24u25) . It turned out that only an asym- 
metric assumption |20j about the mutual screening can correctly describe the 
Onsager resuh [TOIITTIIT^ . 

Due to the widespread usage of T-matrix approximations ranging from 
nuclear physics for both equilibrium p51[?fl[^ and non-equilibrium [^150] 
problems, to the theory of moderately dense gases [3T] and liquid ^He ^32], 
as well as to electron-electron correlations in molecules and solids ^33^34u35l 
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[36ll37] it is desirable to clarify the correct structure of the Bethe-Salpeter 
equation and to specify the conditions under which the corresponding ap- 
proximate forms can be used. 

Here we will present a derivation of modified Bethe-Salpeter equations 
from first-principles which shows indeed that the hierarchical structure of 
correlation functions leads to such asymmetric selfconsistency as necessary 
to describe pairing. We will present the derivation in terms of the causal 
n-particle Green's functions. The cummulant expansion yields a hierarchy of 
correlation functions which leads already to the asymmetric result assumed 
by Kadanoff and Martin in binary correlation approximation. 

We will then proceed and investigate the three-particle correlation and 
will select the processes between two pairs of particles which allow to con- 
sider the repeated collisions. We will show that the systematic cummulant 
expansion leads to the correct subtraction in the channels with condensates 
called singular channels. The idea of derivation follows here the centennial 
overview of many-body approximations by Heinz Puff |38j which is suited for 
nonequilibrium Green's function expansions. Therefore all outlined formal- 
ism holds in nonequilibrium as well as equilibrium. The expansion scheme of 
three-particle correlations will lead to the appearance of two different kinds 
of selfconsistent propagators in each considered channel. These will be the 
channel-dressed and the completely dressed propagators. 

The outline of the paper is as follows. In the next chapter a scheme of 
Green's functions with respect to the cummulant expansion is shortly pre- 
sented. This is performed with the help of the variation technique including 
an external potential in order to ensure the Ward-identity, 0-derivability and 
conservation laws which will be shown in the appendix. Then we consider in 
the third chapter the binary approximation neglecting any three-particle cor- 
relation. This will lead to the Kadanoff- and Martin approximation providing 
the correct gap equation. We discuss then the repeated collisions and how 
it can be avoided by a proper subtraction scheme. Though leading to the 
same gap equation the underlying approximation is beyond the KadanofF- 
Martin form. To justify this subtraction scheme we consider in chapter IV 
the three-particle correlations. We demonstrate that the expansion of three- 
particle correlations in different channels leads to the distinction of two kinds 
of selfconsistent propagators, the channel dressed and the completely dressed 
one. Then we focus on the correlations between two pairs. This provides a 
special contribution from the pairing or condensation channel to the different 
channels contained in the three-particle correlations. It turned out that this 
selfenergy contribution provides exactly the subtracted scheme discussed in 
chapter HI and which was proposed in the literature. In the appendix we 
proof the Ward identity and the «?-derivability in order to ensure conserva- 
tion laws. We summarize and emphasize that the hierarchical dependence of 
correlation care for the avoidance of unphysical repeated collisions itself. As 
a consequence the T-matrix schemes used so far should be revised as soon 
as singular channels appear like in phenomena of condensation or pairing in 
such a way that an asymmetry in the selfconsistency is proposed. 
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2 Decoupling scheme of correlations 

The n-particle causal Green's function 



G(l,2,...,n,l',2',....,n';f/) = 



1 {TS^i...W,,^+...^+) 



(TS) 



(1) 



of a system under the influence of an external field U{1), where numbers 
1,2,... sign cumulative indices like space, time,... coordinates, can be formally 
represented by a generating functional 



G(l, 2, n, 1', 2', n'; U) = d^^, ...d^^, dx, -.dx^QiX. ^] L=^=o (2) 



with 



c;[A,^] = i+;^^ y di...dndv. 



.dn'K...\iG{l...n'-U)r]v...r]n' (3) 



where r;, A are Bose/Fermi-commuting auxiliary fields. 

The Martin-Schwinger hierarchy coupling the one-particle Green's 
function to the two-particle one. 



y2 

^dt, U{1) 



G(ll'; U) = S{1 -l')Ti / dlV{l, 1)G(111'1+; U), 



(4) 



is expressed in terms of this generating functional by 

idt, + ^^u{i)\dx,g = mGT^ I div{i,i)d,,^^^dx,dx,g. (5) 

2m J J 'IT 1 

The upper sign denotes fermions the lower bosons hereafter. It is now useful 
to introduce the correlated n-particle Green's function as the cummulant 
expansion '40| due to a new generating functional Qc, 



g[\,Tf\^ exp 



ith 



n . 

0,[A,r;] =^^y dl...dndl'...dn'Xn..MG,{l...n'- 



(6) 



u)m'-vn'. (7) 



The comparison of ([61), Q and ([3]) with respect to the orders of A, 77 reveals 
the cummulant expansion. The first order reads 



G(ll';[/) = G,(11';{/), 



(8) 
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i.e. the single-particle Green's function equals its correlated part. The sec- 
ond order shows just the separation of the Hartree-Fock term from the two- 
particle Green's function 

G(121'2';C/) =Gc(121'2';t/) 

+ G(ll'; U)G{22'; U) T G(12'; C/)G(21'; U). (9) 

The third order is given by all possible exchanges of three one-particle Green's 
functions together with all possible exchanges of the two-particle correlated 
Green's function and a one-particle Green's function. 

Introducing ([6]) into ([Sj and notating the inverse Hartree-Fock Green's 
function with 



Ghf(12; U) = -iV{l2)G{l2-, U) 

5{l-2) (10) 



idt,+^-U{l)±i I d\V{ll)G{ll+-U) 
2m J 

one can invert the differential equation ([S]) into an integral equation 

dx,gc^i{i) (11) 

with 



X(l) = j dlGYiF{ll\U)mTi J dld2GnFill;U)V{l2) 

+ [d^^^dx,g,~G{22+;U)]dx,g, 
T[dr„^dx,gc-G{12+;U)]dx,gc 

d^^^gc[dx,dx^,gc + dx,gcdx-^g,]} . (12) 

Using ([7]) we see that from the definition the correlated Green's functions 
can be represented by functional derivatives of (jlip . 

Gcin';U)^d,^,dx,gc\^^,,^„ 
Ge(i2i'2'; u) - a,,,, a,,^, dx, 9a. gc I ,^^^0 



=T5„^,a^,,5A,X(l)|^^^^o 
Gc (1231'2'3'; [/) = 9^,, d^^, 8,,^, dx, dx, dx, gc 



X=ri=0 



= d,„d,,_,d,^,dx,_dx,Iil)\,^^^, (13) 

and so on. We drop the notation of the explicit dependence on the external 
potential U in the following. The first equation of (IT3| yields for the one- 
particle correlated Green's function 

Gc(ll') = Ghf(11')T« J dl2GHF(ll)F(l,2)Gc(121'2+) (14) 
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and is nothing else but the Martin-Schwinger hierarchy (|4]) written for the 
correlated parts and in integral form. 

The integral equation for the two-particle Green's function is more in- 
volved and reads from the second equation of ([T^ 

Gc(12l'2') = i J dld2GHF{ll)V{12) |TGc(2i22+l'2') 

-H:;(22+) [Gc(121'2') + G(22')G(11')tG(21')G(12')] 
-1-Gc(i22+2')G(21') + Gc(11'22+)G(22') 

TGc(222+2')G(11') t Gc(221'2+)G(i2') | . (15) 

This equation is graphical presented in figure [T] and is exact so far. Since 
the following algebra is somewhat involved it is easier to perform it in terms 
of diagrammatic presentation. Therefore we design a complete selfconsistent 
one-particle Green's function with a thick full arrow, the interaction with a 
broken line and the two-, three-, and four-particle Green's functions with a 
square containing two, three and four legs. The Hartree-Fock Green's function 
is designed as thin arrow. Corresponding upper signs denote fermions and 
lower signs bosons. The numbering of in and outgoing channels are in the 
direction from bottom to top. Lines without an arrow indicates ends of two- 
and three particle Green's functions which are not connected with a one- 
particle Green's function. 






+ i 



+ I 



Fig. 1 The equation for the correlated two-particle Green's function. The thin 
open arrows are Hartree-Fock Green's functions, the thick arrows are full Green's 
functions. The first line represents the ladder approximation and the second and 
third line exchange channels. 
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From the derivation above it is not hard to see that the property 

Ge(12l'2'; U) = tG,(211'2'; U) = TGe(122'l'; U) (16) 

is ensured. In other words the criterion (A) and (B) of Kadanoff and Baym is 
fulfilled and we have a conserving approximation, i.e. the energy, momentum 
and density is conserved [41]. Please note that the corresponding diagrams 
in the second and third line in figure [1] cancel each other for fermions and 
contact interaction, i.e. they become small for short-range interactions. 



3 Two-particle correlations 

3.1 Binary collision approximation 

As a first very drastic approximation we want to consider the binary correla- 
tion approximation which consists in neglecting the correlated three-particle 
Green's function in ([15]), i.e. the last diagram in figure ([T]). Further we ne- 
glect also the second and third line since they vanish for fermions and short- 
range interaction. This results into the known ladder approximation or Bethe- 
Salpeter equation. Amazingly we have now obtained that the intermediate 
propagators of the Bethe-Salpeter equation have to be considered asymmet- 
rically. Indeed, the vertex 

Gc(12l'2') = J dldM'd2'GHF(ll)G(22)r(i2T'2') 

xG(ri')G(2'2')) (17) 
can be represented in symmetrized form 

r(12l'2') = r(12l'2') =F T(122'l') (18) 

by the T-matrix 

r(12l'2') = iV{12)5{l - l')S{2 - 2') 

+ i ydld2T/(12)G(22)GHF(ll)T(12l'2'). (19) 

The asymmetry appears such that a full Green function Gc is combined 
with a Hartree-Fock Green's function. This asymmetry is maintained if we 
go systematically to higher-order approximations neglecting 4-particle corre- 
lated Green's functions as we will see in the next chapter. It is a result of 
the hierarchical structure of the equations of motion. If one compares with 
the linearized parquet approximations [42 ,43 , the Hartree-Fock propagator 
Ghf in figure ([T]) is replaced by the full one Gc which includes nonphysical 
processes. We will see in the next chapter that the result derived here by 
the hierarchy avoids such non-physically multiple scatterings with the same 
channel. Therefore we consider the asymmetric form in figure [1] as superior 
to the parquet approximation. 
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k k k k 

Fig. 2 The Dyson equation in ladder approximation for the symmetric selfcon- 
sistency (above) and the Kadanoff-Martin asymmetric one derived here (below). 



In fact this asymmetry has been recognized as being necessary to obtain 
the gap equation for pairing [li]. It was first observed by Kadanoff and 
Martin [T] and used later on [H[51I^H51 as an ad-hoc approximation seemingly 
violating the symmetry of equations and consequently violating conservation 
laws. This has remained puzzling since a worse approximation leads obviously 
to better results. Recently it turned out that the repeated collisions 9 with 
the same particle are responsible for this artifact. 



3.2 Derivation of Gap equation 

Let us illustrate this by comparing the Dyson equation with symmetric and 
with asymmetric selfconsistency in figure [21 The T- matrix has poles at bound 
states and at the pairing and becomes separable near these poles, T = =fAA 
where the upper sign is for fermions and the lower for bosons. The two 
different Dyson equations for the propagator with respect to the symmetry 
of selfconsistency is illustrated in the next figure [3] The first symmetrical 



+ 



A , A 
k -k 

-A. A 



k ' k * k -k k 

Fig. 3 The propagator equations resulting from figure [2l selfconsistently (above), 
Kadanoff-Martin asymmetric (below). The upper/lower sign stands for Fermi/Bose 
systems. 



9 



selfconsistent Dyson equation leads to the propagator 

= TT (20) 

which shows no pole and no gap equation. Considering the second asymmet- 
rical Kadanoff-Martin approximation of figure [3] we obtain 

G(^,k) = (21) 

which possesses the typical two-pole structure of the BCS gap equation in 
the case of fermions. Therefore the Kadanoff and Martin approximation is 
superior to the symmetric selfconsistent ladder approximation and appears 
as a consequence of the hierarchical dependencies of correlations as derived 
above. 

What has this now to do with avoiding repeated collisions? Repeated 
collisions of two particles in the same state are unphysical since the particles 
move apart from each other after the collision. Therefore we have to ensure 
that due to selfconsistency such collisions with the same state do not appear. 
In fact as can be seen in figure |31 the selfconsisent Dyson equation does 
not ensure that the momentum of the repeated collisions, p, is unequal to 
the momentum of the incoming particle k. If this would be the case the 
particle will scatter with a particle in the same state again. These repeated 
collisions have to be removed from the T-matrix and the correct gap equation 
appears and the condensate can be described without asymmetrical ad-hoc 
assumptions about selfconsistency. The advantage of eliminating only the 
contributions of single channels as proposed in Refs. [5] and [IS] is that the 
formation of pairs and their condensation can be described within the same 
approximation. Explicitly, we split the selfenergy into different channels, S = 




Fig. 4 The iteration of the selfconsistent Dyson equation of figure [5] leading to 
repeated collision with momentum p, q. 



Sj where we assume the condensation or pairing to appear in the channel 
i called the singular channel. Now we define a subtracted propagator 

G^=G- G^U,G (22) 

where the T-matrix of the channel carrying the pole, Tj, is closed by the 
subtracted propagator Si = TiG^^. Writing the Dyson equation explicitly 
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leads to the full propagator in momentum-energy representation p = (i^, p) 

m 

^ _^ + e + Til 

(t^ + e + Ai)(c^-e-^ii)T^?2 

where Enip) = S(p) — Si(p) and Si2(p) = ^(p) and the upper/lower sign 
stands for fermions/bosons. The expression is nothing else but the Be- 
liaev form for bosons [3^ or the Nambu-Gorkov form for fermions |37] and 
can be written also in matrix notation with off-diagonal elements [12j. 

Though we have shown above that the Kadanoff-Martin approximation 
leads already to a two-pole Green's function (PT|) necessary for the gap equa- 
tion, the alert reader has noticed that the subtraction of unphysical multiple 
scattering events has assumed a different intermediate propagator than the 
Kadanoff-Martin approximation though leading to the same gap propagator 
The Kadanoff-Martin approximation requires to use the Hartree-Fock 
propagator as derived above by approximating the hierarchy of correlations 
at the binary level. For the subtraction of repeated collisions we have used 
instead the propagator written with the help of (1221) 

= Ghf + Ghf('S' — SuF — Ei)Gi^ (24) 

which shows that the subtracted propagator is beyond the Hartree-Fock one. 
How can we derive such a form from the hierarchy of correlations? 



4 Three-particle approximations 

4.1 Generic structure 

For this purpose we go one step further and approximate the equation for 
the three-particle correlation represented by the last line of (fT5|) . It shows the 
coupling to the four-particle Green's function and is quite lengthy. For the 
sake of legibility we abbreviate the interchanges of indices of the correspond- 
ing foregoing expressions by denoting them in the following formula within 
the same kind of brackets, 

Gc(123l'2'3') = T«y"<^lrf2GHF(ll)V(12)|Gc(23212+l'2'3') 

T Gc(221l'2'3')Gc(32+) T (2 o 3) 
+ {Gc(2322+2'3')Gc(ll') T (1 ^ 2) 
TGc(222+3')Gc(3il'2') T (2 o 3) 
-H[TGc(32+)Ge(222'3')Ge(il') + (1 o 2)] t [(2 ^ 3)] 
-hGc(232+l')[Gc(212'3') -I- Gc(22')Gc(13') 
TG,(23')Ge(i2')]} 

T {1' o 2'} + {1' o 2' o 3'} 1 . (25) 
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Now we consider selected approximations by neglecting the four-particle cor- 
relation and selecting special sets of diagrams. First let us show how the 
known channel approximations appear which are the screened ladder, the 
maximally crossed diagrams and the ladder diagrams. Then we will consider 
the pair-pair correlations and their influence on these three channels. 






Go 


+ i 


Go 


+ i 


Go 





Fig . 5 A special set of diagrams for the three-particle correlated Green's function 
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Fig. 6 Definition of the screened potential Vs. 




4.2 Screened ladder approximation 



We choose as partial summation from the first set of diagrams indicated 
in figure[5l The diagrams obtained by interchanging 1' ^ 2' and f ' <-> 2' o 3' 
in figure [S] are contained in (051) as well. Now the iteration of the equation for 
the three-particle correlated Green's function in figure [5] leads to a repeated 
sum in the interaction lines which can be summarized by introducing the 
screened potential of figure |6l This procedure results into the expression 
for the three-particle correlated Green's function as illustrated in figure [71 
Introducing this expression into the last diagram of figure [T] one obtains the 
diagrams of figure [51 We see that the last line renormalizes the single-particle 
propagator if brought to the left-hand side. In such a way the equation of 
the two-particle Green's function in figure [51 can be very much simplified to 
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-G 

c 


p3' — 1 
- 2' = + i . 




+ i 




Gc 







Fig. 7 The set of diagrams of figure [S] when introducing the screened potential of 
figure [6l 



Go 



C3c 



Fig. 8 The diagrams of figure [T] when introducing the three-particle Green's 
function of figure [T] 



which end we introduce the modified propagator 

Gs = Ghf + Ghf{Ss — Sf)Gs 
= Gh + Gh(^.)G, 

= Go + Go{Us + Sii)G, (26) 

where Z's(12) = iVs{12)G{12) is the screened self energy and Z'f(12) = 
iy(12)G(12) the Fock self energy as illustrated in figurelHl We will call in the 



+ I 



^ -I 



Fig. 9 The screened propagator ()26|) . 



following this modified propagator the channel-dressed selfconsistent prop- 
agator to distinguish him from the complete-dressed selfconsistent propa- 
gator. The channel-dressed propagator is determined by the corresponding 
selfenergy understood as the lowest selfconsistent diagram in the correspond- 
ing channel while the complete-dressed propagator includes all higher-order 
crossed terms. The introduction of this channel-dressed (screened) propaga- 
tor results into the final expression of figure [TUl 

Compared to the Kadanoff-Martin approximation which was represented 
by diagrams of figure [1] neglecting the three-particle Green's function, we 
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Fig . 10 The diagrams of figure[8]when introducing the channel-dressed propagator 
(|26|l indicated as thin open arrows. 

now obtain in figure [TU] the same set of diagrams except tfiat the Hartree- 
Fock propagator has to be replaced by the modified one ([26]) and the bare 
interaction has to be replaced by the screened interaction of figure [5] in the 
ladders. Please note that the screened interaction appears with asymmetric 
propagators. 

Again we see that for fermions and short range screened interaction the 
second and third line of diagrams in figure [TU] are canceling mutually and only 
the screened ladder diagram remains. Together with the screened potential 
of figure [6] this is the screened ladder approximation used before [48] except 
that now one of the internal propagators has to be replaced by the screened 
one pUj) . This establishes the asymmetric form derived here as a new result. 

Please note that (j26p looks already like the structure of subtracting un- 
physical repeated collisions ([M)) . however, only the Fock term appears to 
be subtracted from the screened approximation selfenergy in the propagator 
([^S)) since the scheme accounts for it already in the Bethe- Salpeter equation. 



4.3 Maximally crossed diagrams 

A next set of diagrams included in (1251) is summarized by the maximally 
crossed ladders presented in figure [TT] where the diagrams interchanging 1' o 
2' and 1' -o- 2' -o- 3' are contained in ([25]) as well. We proceed with the same 
steps as in the last chapter by introducing the expression of figure [11] into 
the last diagram of figure [T] but defining now the channel-dressed propagator 

Grp = Gq -\- Gq X/Ji Grp 

— Ghf + Ghf(^j' ~ 'S'hf)Gt (27) 
to obtain the diagrams in figure [T^] Here the selfenergy reads 



Sf{n')^T / dl2/'(1121')G(21+) 



(28) 



14 



3' 

2'= -i 
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Fig. 11 The maximally crossed diagrams contained in (|25p. 
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X 





Fig. 12 The diagrams introducing figure [TT] into figure [T] The thin open arrow 
mark now the channel-dressed propagator ((27)l . 



where the symmetrized vertex 

f(1234) = r(1234) =F 7(1243) 
is expressed via the maximally crossed ladders of figure [T31 



(29) 



Fig. 13 The maximally crossed ladder summation. 



4.4 Ladder diagrams 

As a last set we select the ladder diagrams included in (p5)) which are collected 
in figure [TH Also the diagrams interchanging 2 o 3 are contained in (OSl) . 

Introducing the expression of figure [T?] into the last diagram of figure [1] 
with the channel-dressed propagator 

Gt — Gq -t- GqUtGt 

— Ghf + Ghf(^t — Shf)Gt (30) 
we obtain the diagrams in figure 1151 Here the selfenergy reads 

Z't(II') = y"dl2r(112l')G(21+) (31) 

and the symmetrized vertex F was given by (|18p in terms of the ladder 
T- matrix (1191) presented in figure 1161 
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2' = + i — T 

I 

-1' ■• 



G 
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Fig. 14 The ladder diagrams contained in (|25|) . 
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Fig. 15 The diagrams introducin g fig ure 1141 into figure [T] The thin open arrow 
are the channel-dressed propagator H30I) . 



Fig. 16 The ladder summation. 



4.5 Pair-pair correlation 



Now that we have seen how the standard channels appear from the cum- 
mulant expansion of correlations in terms of Green's fmictions with special 
emphasis on the asymmetric propagators we proceed and investigate the pair- 
pair correlation. In fact in (j25p there are diagrams included which describe 
the interaction of two two-particle Green's functions outlined in figure [T71 
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Go 
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Gt 
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1' 

-- 1 
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Gt 



Gt 



Gt 



Fig. 17 The pair-pair diagrams in (|25 



From these diagrams we search only the ones which yield a renormaliza- 
tion of the left-lower Hartree-Fock propagator after iteration in one of the 
above three channels, i.e. the first responsible diagram of figures [5l [Til or [T4l 
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respectively. In such a way we can define a channel-dressed propagator as 
done repeatedly above. It turns out that only the one diagram written in the 
second line of figure [T7] fulfills this task. All other diagrams give repeated 
iterations partially included in the above summations and partially leading 
to new cross diagrams. These diagrams we will not consider here. They have 
been partially considered in Born approximation named as cluster Hartree- 
Fock diagrams [J^ and have been applied to exciton problems [S^ and to 
the first-order superfluid phase transition J5T]. These diagrams describe the 
interaction between the cluster and the single particle, while we concentrate 
here on a genuine cluster-cluster diagram. 





Fig. 18 The selfenergy of the renormalized propagator due to pair-pair correla- 
tions in figure [17] where the kernel K signs Vs — V for the screened channel of figure 
1101 iT for the ma xim ally crossed channel in figure [12] and iT + V for the ladder 
diagrams in figure [TS] The second line is valid only for separable two-particle cor- 
relations, see figure 1191 when using the screened dia grams of figure [8] The third 
line appears using the ladder diagrams of figure 1151 Crosses denotes the inverse 
(amputated) propagator. 



We introduce now this single renormalizing diagram of the last line in 
figure [iTj into the corresponding three-particle ones on the right hand sides 
of figures [5l[TT] or [lU This leads to iterations which sum the interactions and 
result into the channel effective blocks, screened potential, maximally crossed 
vertex or T-matrix vertex. These diagrams are then again introduced into 
the last diagram of figure [Tj as done repeatedly before. In this way we obtain 
an additional renormalization diagram for the channel-dressed propagators 
of (dSJ), (dZ]) and ((Sni) which can be written 

Ga — Gp ~ GpSaGa (32) 

with p = s,T,T denoting the channels. The corresponding selfenergy Sa due 
to the pair-pair interaction is shown in figure [T8| where the kernel K denotes 
the considered channels Vg — V, iT and iT + V correspondingly. 
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Fig. 19 Separation in in- and outgoing channels in the pole of pairing or conden- 
sation. 

We want to consider specially the particle-particle channel represented 
by the T-matrix since in this channel the pairing appears. There we have 



For the singular channel, where the pairing appears, we will see now that 
the desired specific contribution remains. For this purpose we use the fact 
that the two-particle correlation separates near the pairing or condensation 
pole, illustrated in figure [19] and as it was used above in chapter [3T2] Taking 
this into account we see from the screened ladder diagrams in figure [5] that 
besides the Hartree and Fock terms which appear to be subtracted in figure 
I18[ only the third term in figure [5] remains as a connected diagram. Therefore 
we can replace the corresponding form in the second part of figure [18] by a 
full two-particle propagator which results into the second line of figure [TS] 
Using the iteration of the T-matrix channel of figure [T5] once more, we obtain 
the third line in figure [TS] Due to this procedure we obtain the result 



and comparing (|^^ with (l55|) we concluded that Ga — G^. This is exactly the 
subtracted propagator (|24|) proposed before which corrects for the repeated 
collisions with the same state. The derived selfenergy diagram in the first 
line of figure [TS] is therefore one of the main results of this paper since it 
presents the diagrammatic part which leads to the subtraction of unphysical 
repeated collisions in any channel. For the particle-particle channel we have 
obtained the correct subtracted propagator and see that a proper collection 
of pair-pair correlations in the hierarchical expansion of correlations cares for 
the subtraction of such unphysical processes. 



5 Summary and outlook 

We have re-examined the Martin-Schwinger hierarchical structure of correla- 
tions. The expansion scheme of Green's functions has been written with the 
focus on the cummulant expansion form. We have investigated the binary col- 
lision approximation as well as three-particle correlations and have derived 
the corresponding Bethe-Salpeter equations. It appears that the two internal 
propagators of the resulting Bethe-Salpeter equations are asymmetric with 



Ga 




(33) 



(34) 
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respect to self-consistency. Dependent on the considered channel we obtain 
one complete-selfconsistent and one channel-dressed selfconsistent propaga- 
tor. The binary collision approximation leads to one Hartree-Fock and one 
selfconsistent propagator in agreement with the KadanofF and Martin form. 
This leads to the gap equation of pairing while the standard Bethe-Salpeter 
equation with two selfconsistent propagators does not. 

We have proceeded to understand this asymmetry as a subtraction of re- 
peated collision with the same state which is unphysical. Such single-channel 
corrections contribute to diagrammatic sums with a weight of the inverse 
volume and appears therefore only if one has a singular channel carrying 
a condensation with a macroscopic number of occupation. It is shown that 
such a subtraction scheme appears if we consider the correlations of two pairs 
of particles contained in the three-particle correlated Green's function. We 
identify the selfenergy diagram which is responsible for such proper subtrac- 
tion scheme and show that the recently proposed scheme following the idea 
of Soven has its justification in the hierarchical dependence of correlations. 

The here derived selfenergy correction responsible for this subtraction 
scheme is valid in all three considered channels, the screened ladder, the max- 
imally crossed diagrams and the T-matrix channel. Though we had searched 
only for the correction of the latter one in order to obtain pairing and the 
gap equation we suggest that this selfenergy diagram might have an effect 
also to the maximally crossed diagrams. In that channel one describes weak 
localization phenomena. Since it was shown recently ^43, that the parquet 
diagrammatic summation leads not to the expected Anderson localization, 
we suggest that the asymmetric parquet summation which appears due to 
the here derived selfenergy correction might solve these problems with local- 
ization. This investigation is deserved for further work. 
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A Variational technique, Ward identities and ^-derivability 

Due to the coupling of an external potential [7(11') — 5ii/U{l) one expresses the 
two-particle causal Green function G(121'2') = 1/i^ < T<I'i'I'2'I'2'^t > by a varia- 
tion of the one-particl e Green 's function G(12) — 1/i < > with respect to 
the external potential [521153] as 



where the upper sign denotes the Fermi and the lower the Bose functions. Using 
the Dyson equation 



G(121'2') = G(ll')G(22') t 



SGjll') 
5U{2'2) 



(35) 



G"' = G, 



-E-U 



(36) 
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one can calculated the derivative in psp and with the help of the chain rule and 
SG — —GSG^^G, the fluctuation function reads 

L(121'2') = Gc(12l'2') + G(12')G(2l') 
= G(121'2') - G(ll')G(22') 

= TG(12')G(21') t G(13)^P|g(41') 
=TG(12')G(2l')+G(13)|||^i(5262')G(4l'). 

(37) 

Double occurring indices ar e un derstood as integrated over. Equation (|37p is ex- 
pressed graphically in figure [50] 




~ 34 — H 



5G 



56 



'3645 



Fig. 20 Variational representation of the two-particle correlation function (|37|l 



Defining a three-point vertex function 



r(123) = SiiSvi + 



<5i25i3 + ||g|G(46)r(673)G(75) 



(38) 



illustrated in figure 1211 t he two-particle correlation function of figure [20] can be 
represented by figure !^ 

The first line of equation (|38p establishes the Ward ide ntity 54 , 55 , 56 which is 
therefore a co nsequence of the variational expression (|35p and the form of Dyson 
equation ([36]). This Ward identity is also a consequence of the fact that the selfen- 
ergy can be represented as variation of a functional with respect to the propagator 




Fig. 21 Definition of the vertex function (|38p 



20 




1 2 

Fig. 22 Equation for the vertex function. 





Fig. 23 Equation for the vertex function in asymmetric form. 




iX Xr 
+ i 




Fig. 24 Kernel of the vertex function in figure [22] for screened ladder channel of 
figure [TU] 



[39| . Recently it has been proposed [57] to expan d th e variational scheme in terms 
of the kernel represented in the last line of figure 1201 

How does this scheme fit now the asymmetric form presented in this paper? 
We demonstrate this on the example of the screened ladder approximation of figure 
1101 The other channels can be t reated analogously. Closing the upper lines of the 
mirrored diagrams of figure [10] with the channe l-co rrected propagator we obtain 
the same form of the three-point vertex in figure [21] but now translated into figure 



23] with asymmetric propagators. The kernel takes the form represented by figure 
24 which can be verified by e xpa nding the diagrams with the help of Gc of figure 
Tin In comparison with figure [22] the intermediate propagators are asymmetric. 

Using the Dyson equation for the channel-d ress ed propagator with channel 
index i we can replace in one propagator of figure [22] 



G — Gi ~\- GiSiG 



(39) 



and then have to show that 



£■(1234) + =(1235)Gi(56)r,(64) = ^(1234) (40) 
in order to proof F = F. To complete (I40p we use 1 -I- GEi — GG~^ and write 
=■(1234) = H(1235)Gi(56)Gr^(64) 

(41) 



5E{12) _ SE{12) (5G(56) 



<5g(43) (5G(56) 5^(43) 
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where we have introduced the (variational) selfconsistent propagator Q to be defined 
by 



^G(56) 
55(43) 



G,(42)G(61) = G(31)G(52) 



(42) 



illustrated in figure [25l 



5G56 


S-^ 


S^43 





= x 



5 6 1 

Fig. 25 Definition of the variational propagator Q (|4? 



With the help of this variational propagator we can fulfill the relation (|41l) such 
that we have 



r(1234) = 



5E{12) 
55(43) ■ 



(43) 



This, however, has the required form of figure [201 such that (|37l) is fulfilled and 
the Ward identities are completed. Please note that neither G nor d serve as 
prop agators to perform the variation but the variational propagator Q defined by 

Therefore the screened ladder approximation (and the other 2 considered chan- 
nels as well) can be recast into a form of three-point vertex function which obeys 
an integral equation with has the desired form of Ward identity. The difference 
to the traditional treatment is that a full dressed propagator is combined with 
the channel-dressed propagator. The price we pay is that we do have two different 
selfconsistent propagators the channel-dressed one and the complete dressed one. 

Now we can write down the standard ^-functionals for T-matrix, screened lad- 
der and maximally crossed channels but with symmetrical Q propagators. 




Fig. 26 The cf-functional defining the asymmetric screened channel approxima- 
tion. 
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Alternatively we can use H42|) to replace them by the two different selfconsis- 
tent propagators. Consequently, the appearance of two different propagators does 
not violate the #-derivability and the conservation laws. This can be considered 
alternatively by translating the known conserving approximations into forms of 
two different species in the system and interpret in the end one specie as the 
complete-dressed selfconsistent and the other specie as the channel-dressed selfcon- 
sistent propagator. Fo r th e screened ladder channel we declare e.g. a ^-functional 
as illustrated in figure [251 

At the end it has to be remarked that the here proven Ward identity and 
derivability is in a sense a triviality as long as the scheme is kept exact. Any formal 
expansion is 'exact' and will obbey therefore these exact identities. The difference 
between different expansion schemes lays in the different kind of approximations 
'naturally' offered. Therefore at any stage of approximation actually used, one has 
to ensur e th e Ward identity. The introduction of the variational propagator (|42p 
in figure [25] shows that any asymmetric approximation proposed with the scheme 
in this paper can be translated into a symmetric one. However the relation is very 
cumbersome in such a way that any asymmetric form chosen translates into many 
higher order diagrams in the symmetric writing. 
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